Ab initio computation of d—d excitation energies in low-dimensional 
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Using a quantum chemical cluster-in-solid computational scheme, we calculate the local d-d ex- 
citation energies for two strongly correlated Mott insulators, the oxychlorides TiOCl and VOC1. 
TiOCl harbors quasi-one-dimensional spin chains made out of 5= 1/2 Ti 3+ ions while the electronic 
structure of VOC1 displays a more two-dimensional character. We find in both cases that the lowest- 
energy d-d excitations are within the t2 g subshell, starting at 0.34 eV and indicating that orbital 
degeneracies are significantly lifted. In the vanadium oxychloride, spin triplet to singlet excitations 
are calculated to be 1 eV higher in energy. For TiOCl, the computed d-level electronic structure 
and the symmetries of the wavefunctions are in very good agreement with resonant inelastic x-ray 
scattering results and optical absorption data. For VOC1, future resonant inelastic x-ray scattering 
experiments will constitute a direct test of the symmetry and energy of about a dozen of different 
d-d excitations that we predict here. 



CO 



I 

C 

O 

o 



> 

On 
O 



I. INTRODUCTION 



The vast majority of solid-state electronic structure 
calculations is based nowadays on density functional the- 
ory (DFT). For strongly correlated systems, however, 
the wavefunction based ab intio methods certainly pro- 
vide a viable alternative. For this class of materials, of 
which a typical example are the parent compounds of 
the high-T c copper oxide superconductors, it is essential 
to properly describe the multiconfigurational character 
of the many-electron wavefunction. The multiconfigu- 
ration and mult i- reference quantum chemical methods 
constitute here a natural choice. When i) the clusters 
on which such computational techniques are applied are 
large enough and ii) the embedding potential describing 
the crystalline surroundings is judiciously constructed, 
the wavefunction based approach yields results in excel- 
lent agreement with the experimental data. This has 
been recently shown for the sequence of spin-state tran- 
sitions in a d 6 system, LaCo03fi and the <i-orbital elec- 
tronic structure of a number of Cu d 9 oxides^ We used 
in those studies a fully ab initio embedding scheme in 
which the effective embedding potential is constructed on 
the basis of prior periodic Hartree-Fock (HF) calculations 
for the extended crystal, in contrast to other quantum 
chemical investigations based on simpler point charge 
embeddingsi^ Further, as a general recipe, the embed- 
ded cluster that enters the post-HF treatment comprises 
in our approach in addition to the "central" MC»6 octahe- 
dron the nearest-neighbor (NN) octahedra, where M is a 
transition-metal ion. The finite charge distribution of the 
NN octahedra and the tails of "central" Wannier orbitals 
extending to neighboring sites are thus described with 
high accuracy. For the particular case of d-d excitation 
energies in Cu d 9 oxides^ excellent agreement was found 
with state of the art resonant inelastic x-ray scattering 
(RIXS) experiments^ Yet, as compared to earlier quan- 
tum chemical calculations on smaller single-octahedron 



[CuOg] clusters embedded in arrays of point charges, 5 
differences as large as 0.5 eV were observed^ 

In the present work, we verify the performance of our 
cluster-in-solid computational scheme for the case of low- 
dimensional Ti d 1 and V d 2 Mott insulators, in particular, 
the TiOCl and VOC1 compounds. The d-level electronic 
structure of TiOCl has been previously addressed by op- 
tical absorption experiments corroborated with model 
Hamiltonian configuration-interaction (CI) simulations 
of the optical spectra*^ DFT-based investigations^— 
ab initio quantum chemical cluster calculations^ and 
RIXS measurements^ The V oxychloride compound, on 
the other hand, is much less investigated. On the ex- 
perimental side, basic information concerning its d-level 
electronic structure is available from optical absorption 
dataJ^ Here, we perform a detailed quantum chemical 
investigation to better understand the nature and sym- 
metry of the different excited states identified by optical 
experiments 



II. STRUCTURAL AND COMPUTATIONAL 
DETAILS 

TiOCl and VOC1 display the same highly anisotropic 
crystalline structure, with a bilayer network of metal and 
oxygen sites well separated along the crystallographic c 
axis by double CI layers (see Fig. 1). Each metal ion lies 
in the center of a distorted O4CI2 octahedron. Adja- 
cent octahedra share edges along the b axis and corners 
along a. The unit cell is orthorhombic, with space group 
Pmmn (number 59) j 15 i 16 

Our computational investigation implies a sequence of 
several steps. Periodic restricted HF calculations are first 
performed with the CRYSTAL packaged We used the ex- 
perimental lattice parameters and all-electron Gaussian- 
type basis sets. Further post-HF correlation calculations 
are carried out on finite embedded clusters. The same 
type of cluster is employed for both TiOCl and VOC1. 
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It consists of nine MO4CI2 octahedra, i.e., one reference 
octahedron plus eight neighboring octahedra, where M is 
either Ti or V. The reference octahedron in the center of 
this cluster C defines the active region Ca where local d-d 
excitations are explicitly computed. The adjacent eight 
octahedra constitute a buffer region Cb- The role of the 
sites in Cb is to provide support for the tails of Wannier 
orbitals (WO's) in the Ca region. 

The orbital basis for the correlation calculations is a 
set of projected HF Wannier functions. Localized HF 
WO's are obtained with the Wannier-Boys localization 
module= of the CRYSTAL package and subsequently pro- 
jected onto the set of Gaussian basis functions associated 
with the atomic sites of the cluster C. For our choice of 
Cb, the norms of the projected WO's centered within the 
active region Ca are not lower than 99.8% of the original 
HF WO's. Additionally, the HF data is used to generate 
an effective embedding potential for the nine-octahedra 
fragment C. This potential is obtained from the Fock 
operator in the restricted HF calculation and models the 
surroundings of the finite cluster, i.e., the remaining of 
the crystalline lattice. 

While the doubly occupied WO's in the buffer region 
Cb are kept frozen, all valence orbitals centered at ligand 
and metal sites in Ca (and their tails in Cb) are further 
re-optimized in multiconfiguration complete-active-space 
self-consistent-field (CASSCF) calculations. In the lat- 
ter, the ground-state wavefunction and the crystal-field 
excited states at the central Ti/V site are computed in 
state-averaged multiroot calculations. The <i-level split- 
tings are finally obtained by additional single and dou- 
ble CI (MRCI) calculation^ The whole procedure is de- 
scribed in more detail in Ref. |2|. 

The CASSCF and MRCI investigations are carried out 
with the molpro packaged The effective embedding 
potential is added to the one-electron Hamiltonian in 
the CASSCF/MRCI computations with the help of the 
CRYSTAL-molpro interfaced 




TABLE I. Ti 3+ d-level splittings in TiOCl. A non-standard 
local reference system is chosen for the octahedrally coordi- 
nated Ti ion, see Fig. 1, with a rotation of 45° around x. The 
ground-state dy 2 _ z 2 configuration is taken as reference. All 
energies are given in eV. 

Orbital occupation CASSCF/MRCI CASSCF/CASPT2^ 



d 1 2 
y - 

dxy 
dy Z 



0/0 
0.36/0.34 
0.65/0.63 
1.51/1.55 
2.35/2.36 



0/0 
0.29/0.29 
0.66/0.68 
1.59/1.68 
2.30/2.37 



FIG. 1. Layered crystal structure of TiOCl and VOC1. A 
distorted O4CI2 octahedron is highlighted on the right hand 
side. 



a Rcf. 12 



We used O and CI basis sets of triple-zeta quality as 
designed by Towler et al^ and Prencipe et al£2- For the 
Ti 3+ and V 3+ species, we employed the triple-zeta s and 
p basis sets of Mackrodt et al^ and Dovesi et al£^, re- 
spectively, plus triple-zeta d basis functions as designed 
by Ahlrichs et al~H 



III. RESULTS AND DISCUSSION 
A. TiOCl 

The magnetic properties of TiOCl have been the topic 
of extensive investigations over the last decade. The 
largest intersite interactions are antiferromagnetic (AF), 
occur between NN Ti 3+ 5=1/2 sites along the b axis, 
and give rise to AF spin chains along that direction. The 
interchain couplings are also sizable and responsible for 
geometrical frustration. 

The doubling of the unit cell in the low-temperature 
nonmagnetic phas o 27 ' 28 suggests a spin-Peierls dimerized 
ground state at low temperatures. The observed devi- 
ations from conventional spin-Peierls behavior, e.g., the 
existence of two distinct phase transitions at 91 and 67 K, 
were attributed to the above mentioned frustration of the 
interchain interactions 4- 

As concerns the electronic structure of TiOCl, DFT 
calculations within the local density approximation 
(LDA) indicate sizable occupation numbers for each of 
the Ti t2 g levels^ - — Early interpretations of the unusual 
physical properties of this material thus invoked the pres- 
ence of strong orbital fluctuations.— Inclusion of local 
Mott-Hubbard type correlations within dynamical mean- 
field theory (DMFT) and LDA+DMFT calculations ef- 
fectively pulls one of the t?, g levels below the other two 
components and reduces the occupation of the latter to 
nearly O i 10 ' 11 On the theoretical side, this is confirmed by 
earlier quantum chemical cluster calculations^ as well by 
our present results. 

In Table I, CASSCF and MRCI results for TiOCl are 
listed. The active orbital set contains all five 3d functions 
at the central Ti site and the singly occupied Zd y i_ z i or- 
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bitals at the eight Ti NN sites. We considered high-spin 
intcrsitc couplings, i.e., a ferromagnetic (FM) configura- 
tion of Ti d spins. Although the WO's at the atomic 
sites of C are derived from periodic restricted HF calcu- 
lations for the Ti 3d 1 electron configuration, the embed- 
ding potential is obtained by replacing the Ti 2+ 3d 1 ions 
by closed-shell V 2+ 3d 2 species. This is a good approxi- 
mation for the farther 3d-metal sites, as the comparison 
between our results and RIXS data shows. An exten- 
sion of our embedding scheme toward the construction 
of open-shell embeddings is planned for the near future. 
In addition to the active Ti 3d functions in Ca and Cb, 
all O 2p and CI 3p orbitals at the central octahedron are 
included in the CISD calculations. 

As expected for a d 1 system, the MRCI treatment 
brings minor corrections to the d-d CASSCF excitation 
energies. Our results are compared in Table I with earlier 
CASSCF and CASPT2 (complete- active-space second- 
order perturbation theory) data from Ref. Il2i . A major 
difference between the present and earlier quantum chem- 
ical calculations is the size of the embedded clusters, i.e., 
nine octahedra here and one octahedron in Ref. Il2i . In 
contrast to other systems, e.g., layered Cu d 9 oxides such 
as La2CuC>4 and S^CuC^Chi where differences of about 
0.5 eV were found as function of the cluster size (see the 
discussion in Ref. d) , the agreement between the two sets 
of results is good for TiOCl. 

Good agreement is also found with RIXS and optical 
absorption experimental data, with deviations not larger 
than 0.1 eV. The only exception is the d y 2_ z 2 to d x 2 ex- 
citation, which is predicted to occur at 2.36 eV in the 
quantum chemical calculations and observed at about 
2 eV in RIXS, see Table II. In the optical absorption 
measurements this particular transition as well as the 
d y 2_ z 2 to d xy excitation were not identified, which is re- 
lated to the different selection rules for RIXS and optical 
absorption. 

We recall at this point that in a first approximation 
the total energy of a given state within the d n manifold 
is a sum of a crystal- field contribution, i.e., an on-site 
crystal-field splitting, and a magnetic term (see also the 
discussion in Refs. and0). For the ground-state con- 



TABLE II. Ti d-d excitation energies in TiOCl, quantum 
chemical vs. experimental results. A J(, In 2 term was here 
substracted from each of the RIXS values reported in Ref. H3 . 
see text. All energies are given in eV. 

Orbital occupation 

^xy 
&xz 

(J y z 

d^2 

a Ref.QJ. 
b Ref.Q. 



MRCI 


RIXS^, 


Optics^ 











0.34 


0.32 




0.63 


0.55 


0.65 


1.55 


1.41 


1.50 


2.36 


2.01 





figuration of TiOCl, the NN spin interactions along Ti 
chains parallel to the b axis are quite large, J& « 57 meV,— 
and AF. From the exact Bethe-ansatz solution of the one- 
dimensional Heisenberg Hamiltonian, 30 the AF ground- 
state stabilization energy is J In 2. On the other hand, 
from overlap considerations, the (super)exchange with 
the NN Ti d y 2_ z 2 spins is either zero or much weaker for 
the crystal-field excited states. Since the quantum chem- 
ical calculations were performed for a FM cluster^ for 
a meaningful comparison between the MRCI and RIXS 
data, we subtracted in Table II from the relative RIXS 
energies reported in Ref. [H a term In 2 rs 0.04 eV rep- 
resenting the energy stabilization of the AF ground-state 
with respect to the crystal-field excited states. 



B. VOC1 

In the family of correlated transition-metal oxide ma- 
terials, the vanadium compounds display an impressive 
variety of vanadium valence states, crystalline structures, 
and physical properties. The valence configuration, for 
example, may vary from V 2+ d 3 in VOi-^ 32 - 33 to V 4+ d 1 
in CaV 2 0!j2i and NaV 2 5 ^ While the latter vanadates 
display a layered lattice with ladders of O-bridged V sites, 
VOi- x has a rocksalt crystalline structure. At the level 
of ab initio wavefunction calculations, the V oxides are 
rather unexplored. To our knowledge, only the ladder 
vanadates have been investigated by advanced quantum 
calculations 

Besides having two electrons instead of one within 
the t2 g set of orbitals, VOC1 also displays in compari- 
son to TiOCl a more pronounced two-dimensional char- 
acter as concerns the electronic structure and magnetic 
properties.— In Table III, V d-d excitation energies for 
the VOC1 compound are given. Two different sets of 
CASSCF calculations were carried out, for two differ- 
ent electron configurations at each of the V NN's. Pre- 
liminary CASSCF calculations were first performed with 
eighteen orbitals in the active space, i.e., two t% g orbitals 
at each of the V sites in the nine-octahedra cluster. As for 
TiOCl, we considered high-spin intersite couplings, i.e., 
a FM cluster. It turns out that the dominant ground- 
state configuration at a given V site is d\_ z2 d xy . The 
crystal-field excited states for the central V ion are then 
obtained in state-averaged multiroot computations where 
a group of three additional orbitals was added to the ac- 
tive space. To be sure that the three additional orbitals 
are localized at the central V site and to allow only for 
on-site d-d excitations, the occupation of the NN d y 2_ z 2 
and d xy functions was restricted to 1. Crystal-field split- 
tings obtained from such SCF calculations are listed in 
the third column of Table III. To avoid complications re- 
lated to states that display low-spin intersite couplings, 
only the high-spin on-site excitations were computed in 
this case. 

We note that the ground-state electron configuration 

in the periodic restricted HF calculation is d , 2 3 <iL, but 

y — * -Ly 
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TABLE III. V d-d excitation energies in VOC1. A non-standard local reference system is chosen, see Fig. 1, with a rotation of 
45° around x. The dominant ground-state configuration is d y2 _ z2 d xy , with a weight of 81%. For each of the excited states, the 
weight of the dominant configuration is not smaller than 64%. All energies are given in eV. 

Dominant configuration CASSCF^, CASSCF^ MRCI^ Optics^ 



^2gi S—l 


d 1 y 2_ z id^ y 

dxydxz 




0.36 
0.43 


t\g e gi S=l 


dxydyz 

d 2_ id i 

A ,1 

^xz yz 

d y 2_z2dyz 

d X zd x 2 

d X yd x 2 


1.61 
1.74 
1.75 
2.85 
3.53 
3.53 


tig, S = 


d 2 y2- z 2 ^ 
dy2_ z 2dxy 
d 2_ z 2d xz 

^xy^xz 














0.36 


0.34 


0.3 


0.44 


0.46 


0.4 


1.60 


1.73 


1.5-1.9 


1.71 


1.82 


1.5-1.9 


1.73 


1.88 


1.5-1.9 


2.82 


2.55 




3.52 


3.17 




3.53 


3.19 




1.27 


1.47 


1.1-1.4 


1.28 


1.47 


1.1-1.4 


1.56 


1.75 


1.5-1.9 


1.79 


1.97 


1.5-1.9 


1.81 


2.00 


1.5-1.9 



a High-spin V d 1 2 _ z2 d xy NN's, see text. 
b Closed-shell V d 2 2 „ NN's. 
c Ref.Qj. 



in the construction of the HF embedding potential we 
imposed a closed-shell V d^ 2 _ z2 configuration for the re- 
maining part of the lattice. In a second set of CASSCF 
calculations, we imposed a closed-shell d^ 2 _ z2 electron 
configuration for the V NN's in the Cb region of the clus- 
ter as well. Minor variations were found as compared 
to the CASSCF d-d excitation energies with high-spin 
d} 2 ,od\,, NN's, i.e., less than 0.01 eV for the lowest 

y — z -i* y 

three roots and less than 0.03 eV for the higher roots. 
The CASSCF V crystal-field splittings for closed-shell V 
NN's are given in the fourth column of Table III. The lat- 
ter CASSCF wavefunctions were used as reference wave- 
functions in further MRCI calculations. Those MRCI re- 
sults are listed in the fifth column of Table III. Only the 
V 3d, O 2p, and CI 3p orbitals at the central octahedron 
were correlated in the MRCI treatment. 

The agreement between our quantum chemical data 
and optical absorption results reported in Ref. [HI is very 
good, especially for the lowest two excitations (see Table 
III). The small difference between those excitation ener- 
gies indicates that the splitting between the d y 2_ z 2 and 
d xy levels, singly occupied in the ground-state configu- 
ration, is also small, about 0.1 eV. On the other hand, 
the energy separation between the d xy and d xz levels is 
substantially larger, 0.3 eV, and indicates that orbital 
fluctuations are unlikely to be important in VOC1. As 
concerns the S — singlet excitations, they acquire very 
low intensity in the optical spectra. On the experimen- 
tal side, a better characterization of those states call for 
high-resolution RIXS measurements. 



IV. CONCLUSIONS 

With a quantum chemical cluster-in-solid computa- 
tional scheme, we have determined the d-d excitation 
energies of TiOCl and VOC1. The ground state configu- 
ration is t\ g for Ti and high-spin t\ g for V. The lowest- 
energy d-d excitations are for both materials within the 
tig subshell, starting at 0.34 eV. We conclude that there- 
fore orbital degeneracies are lifted to a similar and signif- 
icant extent in the two systems, which excludes the pres- 
ence of strong orbital fluctuations in the ground state. In 
the vanadium oxychloride, spin triplet to singlet excita- 
tions start at an excitation energy of 1.47 eV and above. 
The computed <i-level electronic structure and the sym- 
metries of the wavefunctions are in very good agreement 
with RIXS results and optical absorption data for Ti- 
OCl. For VOC1, future RIXS experiments will constitute 
a direct and stringent test of the symmetry and energy of 
about a dozen of different d-d excitations that we predict 
here. 
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